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1 Introduction 



It is commonly believed that QCD at high temperature T ^> Aqcd is in a deconfmed phase, 
that is the system is composed by a gas of quarks and gluons interacting with a coupling 
oqcd(T) <C 1- Therefore perturbation theory is reliable at energy scales of order T. This 
picture runs into problems when one probes the softer scale gT, where contributions of the 
same order are generated at every order in perturbation theory JTJ. At high temperature 
the relevant gT contributions are obtained by resumming the one-loop Feynman diagrams 
(Hard Thermal Loop resummation 0). 

However, in non-Abelian gauge theories there is another major obstacle in performing 
perturbative computations. Since the static magnetic fields are not screened at HTL level, 
one finds wild infrared divergences at order g 2 T, scale in which one expects such a magnetic 
mass to be dynamically generated by the interacting theory. Therefore there is a need for 
new resummation methods. 

In ref. [[J we proposed a possible approach. The basic idea is to consider a "cold" T = 
quantum system and add thermal fluctuations at lower and lower frequency scales, until the 
non-perturbative scale g 2 T is reached. In this way the system is thermalized at larger and 
larger distances until the scale g 2 T is reached "smoothly" . This corresponds to performing 
a coarse-graining a la Wilson |Q for the thermal fluctuations. We will call such a method 
Thermal Renormalization Group (TRG). 

This approach is inspired by continuum Wilson or "Exact" Renormalization Group 
(ERG) in which one starts with the microscopic ultraviolet action and adds quantum 

modes with frequency larger than a coarse-graining scale A. By continuously varying A one 
gets the ERG flow equations in the form of partial differential equations in A, which can 
be integrated to obtain the full quantum theory at A = 0. The main difficulty associated 
to the application of ERG to T = gauge theories is that it is not possible to introduce 
the infra-red cut-off A in a gauge-invariant way. Therefore the Slavnov- Taylor identities of 
the theory are violated and one can only hope to recover them once the cut-off is removed 

One encounters similar difficulties when the ERG is applied to thermal field theories. 
If one tries to perform a coarse-graining of both quantum and thermal fluctuations, the 
ST identities are necessarily violated by the coarse-graining scale A and it is much more 
difficult to extract gauge-invariant quantities (see for instance [|i4|). 

However, as explained above, the problem is simpler if one is interested in resumming a la 
Wilson the thermal modes only, namely in TRG. One assumes that the T = renormalized 
quantum field theory is known (from experiments, perturbation theory, lattice or other 
approximation methods) and adds thermal fluctuations for frequencies \k\ > A. In this 
way one interpolates between the renormalized theory at T = (A = oo, no thermal 
modes have been integrated out) and the theory in thermal equilibrium at temperature T 
(A = 0, all the thermal modes have been integrated out). 

For a finite value of the cut-off A one is describing a physical, i.e. renormalized, system 
out of thermal equilibrium, with a A-dependent density matrix p(A). From a physical point 
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of view, the better status of gauge symmetry in this approach can be easily understood. 
Indeed, the expectation values and the Green functions in an out-of-equilibrium system 
differ from those at T = (and from those in thermal equilibrium) because of the different 
weights of the physical states in the thermal average. But, as long as no unphysical states 
are introduced, the BRS variation and the thermal averaging are commuting operations and 
the Slavnov- Taylor identities are completely analogous to the ones at T = 0, provided the 
vacuum expectation values are interpreted as thermal averages with respect to the density 
matrix p(A) (if 



From a technical point of view, varying the density matrix corresponds to changing the 
temporal boundary condition that has to be given in order to invert the quadratic part in 
the Lagrangian and define the free propagator. The Lagrangian itself is not modified and 
retains all its symmetries. This is to be contrasted with the usual formulation of the Wilson 
RG in which the introduction of the cut-off breaks the BRS invariance of the Lagrangian 
explicitly @, 0-0, [14 . 



In order to implement this modification of the temporal boundary conditions we need 
a real-time formalism suited to describe a system with an arbitrary initial density matrix. 
Therefore we find it appropriate to use the Closed Time Path (CTP) method [l5 [| — [l8 p . 



which is a general formalism to deal with out of equilibrium problems. The CTP technique 
is a generalization of the Real Time (RT) ]19[ formulation of equilibrium thermal field 
theories and involves a doubling of the degrees of freedom. 

In CTP the information about the initial state is encoded in source terms concentrated 
at t = to — * — oo. By invoking an adiabatic switching-off for the couplings, such information 
is contained as boundary conditions in the free propagators. By giving equilibrium bound- 
ary conditions, one recovers the path-integral representation for the generating functional 



of the usual Real Time equilibrium thermal field theory given, for instance, in ref. [24 
We modify the boundary conditions by introducing a momentum-dependent temperature 
/3(|&|) such that the path-integral representation for our coarse-grained generating function- 
al can be obtained by substituting, in the equilibrium path-integral, the Bose distribution 
for a cut-off one 

N b (\k \) - N b (\k \,A) = - A) N b (\k \) . (1) 

By varying with respect to A we will obtain the TRG flow equations for the generating 
functional and for all the Green functions of the theory. 

Integrating the flow equations from A = oo [] to A = 0, the resummed Green functions 
to all orders in thermal perturbation theory are obtained. 

Another advantage of using the TRG in studies at high temperature is that the boundary 
condition for the flow equations at A ^> T is just the full quantum theory at T = 0, 
which is directly related to the physical parameters. The presence of the infra-red cut-off 
allows a reliable computation of the initial conditions in the framework of zero temperature 
perturbation theory. 

A good and clear control over the initial conditions is indeed a necessary requirement 



^dn practice, due to the Boltzmann suppression of high-frequency modes, an initial A » T will be 
required. 
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when one is interested in small thermal effects such as the generation of a magnetic mass 
for the gluon, which we expect to be 0(g 2 T). On the other hand, using as initial condition 
the bare theory at some ultraviolet cut-off, as one is forced to do in the ERG, the connec- 
tion between the physical parameters and the initial ones is obtained only after the pure 
quantum corrections have been included via some approximation of the flow equations at 
T = 0. 

There is another advantage in using the real time TRG instead of ERG at finite tem- 
perature in the imaginary time. In a non-interacting thermal gauge field theory it is not 
necessary to give thermal fluctuations to the unphysical degrees of freedom, such as ghosts 
and longitudinal gluons (only physical transverse particles are in contact with the thermal 
bath and get heated). Therefore the Feynman rules to construct perturbation theory can 
be simplified in such a way that only the propagator of transverse gluons receives ther- 
mal boundary conditions. This is true both in physical gauges (such as axial pOfl-p^l) 



and covariant gauges [23[. In this way only the transverse modes of the gluons have to 



be coarse-grained and one gets simplified TRG flow equations. This approach is denoted 
as TRG in physical phase space and is very useful if one wants to exploit TRG as a real 
computational tool. 

The paper is organized as follows. In Sect. 2 the CTP formulation for an SU(iV) gauge 
theory is briefly discussed. In Sect. 3 the TRG flow equations are derived and the kernel is 
discussed. In Sect. 4 we derive the ST identities and show that they are the same as those 
for the theory without the cut-off on the thermal sector. In Sect. 5 the computation of the 
electric and magnetic screening masses at order g 2 T in the TRG formalism is discussed. 
The appendices contain some technical details. 



2 Closed Time Path for pure gauge SU(iV) theory 

In this paper we will specialize to the case of an SU(iV) Yang-Mills theory. The inclusion 
of fermions and other matter fields poses no conceptual difficulties. 

The action is 

S = J d A x{--^F^-F^)-±-F 2 --c-M- c y (2) 

where 

F; v (x) = d,A a u - d v Al + g(A,A A v f , (A, A A v f = e ahc A\A C V , 

■ F„ v = F* u F; u , D„c = d„c + gA )1 Ac; 

^[A^} is the gauge-fixing term and Ai the Fadeev-Popov operator. The metric tensor 
(1,-1,— 1,-1) and the colour indices will be understood in the following. 

The action (|2|) is invariant under the BRS transformations 

5A^ t = r]D ll c, Sc=-7]^cAc, 5c = -—^[A^ , (3) 
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with 7] a Grassmann parameter. In many formulae the fields and corresponding sources 
will be denoted by 

$ = (A, c, c) , J = (j, x, -X) • (4) 

As we stated in the introduction, our formulation of the TRG is based on a deformation of 
the density matrix of the system implemented by introducing an infra-red scale A such that 
for A — ► oo the system is at zero temperature and at A — ► it is in thermal equilibrium 
at a temperature T. For intermediate values of A, only the high-momentum modes with 
\k\ > A are in thermal equilibrium, whereas the low- momentum ones are frozen at T = 0. 

In order to describe this system for a generic value of A, we shall use the CTP formalism 
introduced in [|16|| , which is designed for a system prepared at some initial time to with a 
generic density matrix p. In our applications the limit t — > — oo is always understood. 

The partition function can be represented as a path integral in which the time argument 
of the fields takes values on a path going from an initial time to to t = +oo along the real 
time axis, and then back to to — infinitesimally below it. Indicating with $ x and $ 2 the 
fields with time arguments on the upper and lower pieces of the time path, respectively, the 
partition function is given by (for the proof of this and more details on the CTP method 
the reader is referred to 



, d 4 p 
Tr J ■ $ ' 



Z[J\ J 2 , p] = J ($\ t = t |p|$ 2 , t = t ) exp* {§ Tr $ ■ D, 1 ■ $ + 5 int [$] + Tr J ■ $} , 

(5) 

where 

i Tr $ ■ Dq 1 ■ $ = £ f^- {§4(p) [G ^(p)] iJ Ai(-P) + c>) [So-^F c*(-p)} , 

^^(-F^W + ^l-^W+^-rix^p) - 1-2} 
and Si n t[$] is the bare interaction action 

5 int [$]=5 int [<l> 1 ]-5int[$ 2 ]. (6) 

Go M ^ and So are the gauge and ghost free propagator, respectively (the 2x2 CTP matrices 
are indicated in bold face). 

In order to define the BRS transformation at the quantum level, source terms for the 
composite operators appearing in @, that is 

wJUV) 1 - f^Ac 1 - 1-2, (7) 

have to be added to (|5|). In the following the dependence on the sources and v will be 
understood. 

In order to construct a perturbation theory for the path integral @, one assumes that 
initial correlations at t = to have a small effect on the system at much larger times. This 
hypothesis is equivalent to adiabatically switching off the interactions for t — > to [21|. In 
this way the initial states |<3>\ t = to) (i = 1, 2) are the physical states of the non-interacting 
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theory (this requires special care in the case of a gauge theory, where the Hilbert space is 
larger than the physical space — we will come back to this issue in the following). Moreover 
the initial density matrix describes a free system and can then be taken diagonal, that is 



p~exp l-jd 3 k(3(k) J2 a \ a k > ( 8 ) 

phys. states J 

where as usual is the annihilation operators for the free initial state \<&(k),t = t ) of 
momentum k (for simplicity we choose the same "temperature" f3(k) for all modes). The 
density matrix (|8p is a generalization of the equilibrium density matrix, which is obtained 
for (3{k) = /3(J,ft. The matrix elements of (|§D at t = t can be represented in the following 
way g7|: 

(&,t = t \p\<S>\t = t ) = exp{z± J K ij &&} . (9) 

Therefore the statistical information on the initial state is all contained in the non-local 
boundary source K^(to,x;to,x'), which is concentrated around t = to. is just the 
boundary condition in time needed in order to invert the d'Alembertian operator and to 
define the free propagator. Notice that with a diagonal density matrix, as in (§), the 
boundary sources K % i can be computed exactly (see [HJ and below). 



Our goal is a path integral representation of a system with modes thermalized for scales 
\k\ > A. This can be obtained, according to ([[]), by choosing f3{k) = (3 A {k) 1 where 

!/3a>|£| for \k\ > A 
(10) 
+00 for \k\ < A 

In this way the boundary source K l i becomes A-dependent and the partition function of 
such a system can be written as 

Z A [J\ J 2 , p] = J exp i {i Tr $ • D^ A • $ + S int [$] + Tr J ■ $} , (11) 

where the appropriate source term K A has been absorbed in 

DoX = Do 1 + K A . (12) 

In general, the free propagator (|12"D can be computed by performing the same manipulations 
as are done in the well-known equilibrium case /3(k) = constant [|I7|]. Therefore in a scalar 



theory 0] D\ is obtained from the one in the real time equilibrium formalism of Niemi and 



Semenoff with the obvious substitution (3 — > (3 A (k). This is equivalent to replacing the 



Bose-Einstein distribution function with the cut-off one as in (JTJ) . 

When dealing with a gauge model, there is a further complication, namely, after gauge 
fixing, the Hilbert space and the set of physical states no longer coincide. As a consequence, 



2 One could also generalize the cut-off Bose distribution by using a smooth cut-off function 0(|fc|,A) 
(0(|fc|, A) — ► 1 for \k\ ^> A and 9(|fc|, A) — > for \k\ <C A). In the practical applications, however, we will 
use a Heavyside step function. 
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the partition function cannot be generally defined as the trace over the whole Hilbert space 
of the density matrix p, but is instead given by 

Z{/3) = Tr [Pp] , (13) 

where P is the projector onto the physical subspace. 

We now consider separately the case of covariant and axial gauges. 



2.1 Covariant gauges 



In a covariant gauge, Fadeev-Popov (FP) ghosts and unphysical degrees of freedom populate 
the Hilbert space, and the projector P is non-trivial. This problem is of course closely 
related to the question of which temporal boundary conditions should be given to the FP 
ghosts and to the unphysical degrees of freedom of the gauge bosons, which could never 
come into equilibrium with the thermal bath. 

In equilibrium thermal field theory there are mainly two approaches to handle this 
problem. The standard one (which we call the Bernard-Hata-Kugo approach) consists in 
giving periodic boundary conditions to the unphysical degrees of freedom [^, [2B| . The other 
one (which we call the Landshoff-Rebhan approach) consists in giving thermal boundary 
conditions only to the physical (transverse) degrees of freedom |23|] . We now describe how 
to generalize them to our TRG case. 

(i) Bernard-Hata-Kugo approach. The basic observation behind this approach is that 
the partition function can be rewritten as 



Z{/3) = Tr 



(14) 



where Q c is the FP ghost charge. Moreover, the thermal average of any gauge-invariant 
operator O, 

(O) =Tr[Pp(9]/Tr[Pp] (15) 

is equal to 



(0) = Tr [e- nQc pO\ / Tr 
where Qbrs is the BRS charge. 



if [O,Qbrs] = 



(16) 



Since any physical observable must be BRS-invariant, Hata and Kugo proposed to 
modify the definition of a thermal average from eq. (15) to fll6[). This gives of course Green 
functions different from those obtained by using (|T5j), but the result for physical quantities 
is unchanged. With this definition one can now derive path integral representations for the 
partition function and the Green functions by replacing p with pexp(— nQ c ) and enlarging 
the set of states |<3?\£o) in © to the full Hilbert space. Strictly speaking, Hata and Kugo 
considered the case of an equilibrium density matrix p. However, as is shown in appendix 
A, it is straightforward to generalize their proof to the case in which p has the form 



p ~ exp 



d 3 kf3{k) 



E(A)t (A) . i-t i t- 

4 4 + 2 c h - W% 

A=0 



(17) 
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where and are the annihilation operators for the free state of momentum k for 

gauge bosons with polarization index A, ghost and antighost, respectively. In the special 
case of thermal equilibrium, (3{k) = /3u,^ this corresponds to assigning periodic boundary 
conditions to all the bosonic degrees of freedom and also to the FP ghosts.0 This is the 
approach that has been traditionally followed in the applications of perturbative thermal 
field theory to gauge theories. 

As explained in the previous section the TRG formulation is obtained by choosing the 
"temperature" j3{k) as in (TO). The computation of the boundary source K\ is completely 
analogous to the equilibrium case. The free TRG gauge propagator is then the usual real 
time equilibrium gauge propagator with the Bose distribution replaced by the cut-off 
one in eq. ([!]), that is 



G 0Aiia/ (k) = - (A- + B ^ + a^V) (Q[Ao(*0] + [ A o(*0 ~ A*(A;)] N b (\k \,A) B 
where we introduced the notations 

f(k) (f(k) — f*(k))8(—k 



Q[/(*0] 



(f(k)-f*(k))6(k ) 
1 1 



-/*(*) 



B 



(19) 



(20) 



and the free Feynman propagator 



Ao 



1 1 



1 

k 2 + is 



The usual basis of tensors 
A^y (j^y — 

B 



1 



(ku) 2 — k 
-k 2 I 



2 [(ku)(Uftk v + Uukfj) - k^ky - k 2 u^u v } 



(ku) 2 -k 2 \ 11 



(ku)kA I (ku)k h 

U„ rr^ - 



k 2 



k 2 



(21) 



a 



flV 



D 



ft]' 



2[(ku) 2 - k 2 } 
kfi k v 



(ku)k p 



u v - 



(ku)k h 



satisfy the multiplication properties 



A- B = A-C = A-D = , A + B + D = l, 
A - A = A , B ■ B = B , C-C=\(B + D), 
TtC-B = TrC ■ D = . 



(22) 



3 This last prescription, formulated for the first time by Bernard P5[ , comes from the fact that expression 
( |l4| ) can be regarded as the partition function of a grand-canonical ensemble with purely imaginary chemical 
potential ^fp = wr//3 for FP ghosts 
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In the thermal reference frame we have = (1,0,0,0), and therefore the only non- 
vanishing components of A are the spatial ones, giving the transverse tensor 

kik ■ 

Ay = —5ij + -=?— , (23) 

\K\ 

In the £ — > limit, we have 

Ao(Jfe) - A*(k) — > -2t7c5{k 2 ) , (24) 

that is, the thermal part of the tree-level propagator involves on-shell degrees of freedom 
only (only real particles belong to the thermal bath). The free TRG ghost propagator is 
the same as for a scalar massless particle (notice again the cut-off Bose distribution): 

Soa(*0 = -Q[A (fc)] - [Ao(fc) - Ao(fc)] N b (\k \,A) B , (25) 



(ii) Landshoff-Rebhan approach. A second possible formulation of thermal gauge theo- 
ries in covariant gauges is obtained by restricting the space of initial states |$\£o) to the 
physical one, as required by the original definition of the partition function, eq. flT3|). One 
has to recall that the initial state is given by free particles. In the case of an SU(iV) gauge 
theory, this means that only the gauge boson states with spatially transverse polarization 
will contribute to the trace, and consequently only the transverse components of the gauge 
fields will couple to .fT-sources. This corresponds to choosing an initial density matrix p 

P ~ exp | - J d 3 k (3(k) £ 4 A)t 4 A) | , (26) 

where the sum over A now runs over transverse polarizations only. 

In the case of thermal equilibrium, this choice leads to the Landshoff-Rebhan formula- 
tion of Real Time thermal field theory [^3] . 

The TRG formulation is again obtained by using the cut-off Bose distribution. The 
TRG free propagator is then 

G 0Atia/ (k) = -A^ P A [A (fc)] - (B^ + ar£V) Q[A (fc)] , (27) 

where we introduced the compact notation 

P A [A ] = Q[A (fc)] + [A (fc) - A* {k)]N b (\k \,A)B . 



The free ghost propagator is simply the T = one: 

S„(A0 = -Q[A (A;)] . (28) 

Compared to the standard TRG propagators in fll8] ) and (p5|), we see that only the trans- 
verse component of the free gauge propagator proportional to A^ v gets a thermal structure, 
which incidentally is gauge-parameter- independent. 

In the practical applications of the TRG in covariant gauges, the formulation a la 
Landshoff-Rebhan turns out to be more convenient, due to the absence of thermal ghost 
loops. 
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2.2 Temporal Axial Gauge 



In physical gauges, such as axial or Coulomb ones, the Hilbert space contains no unphysical 
states, and then P=l. In particular the ghosts decouple. 

The derivation of a path integral representation in these gauges then goes on analogously 
to what is done for the Landshoff-Rebhan case in covariant gauges, with only the physical 
transverse gluon modes included in the density matrix (|8]). 

Here we are mainly interested in the Temporal Axial Gauge (TAG), which is obtained 
by taking T a = u^A ^ in (|j), with u M = (1, 0, 0, 0), and letting a —* oo. 

The TAG is not a "good" gauge, in the sense that the gauge is not completely fixed by 
the choice Aq = 0. This gauge freedom manifests itself in the spurious l/k poles in the 
longitudinal part of the free propagator. In the framework of perturbation theory, different 
prescriptions have been proposed to circumvent these singularities. In particular in the real 
time formalism at T ^ the issue has been studied by James in ref. pilfl . As we will discuss 
in subsect. 4.2, when computing the gauge-independent pole masses using the TRG, all the 
terms containing these spurious poles cancel. So in our applications it is not necessary to 
fix a particular prescription. 



The free TAG TRG propagator in covariant form is given by |2D 

Goa,^ = -A- Pa[Ao(A;)] - Q' [H^A (k)] , (29) 

where 

H^ u = B^ u — C M „ , (30) 

and Q' and k' Q has been chosen for definiteness according to James's prescription, i.e. 

( f(k) (f( k )-f*(k))e(-k 3 ; 

wm = 

V (f(k) - r(k))d(h) -nk) 

and k' = ko + ie/k^. 



(31) 



3 TRG evolution equations 



By taking the derivative with respect to A of ( JTT| ) we obtain the evolution equation for 

z A [J], 

A§- A Z M = - iW ^{ |A ^.^l} (32, 

In deriving the second equation the ghost propagator has been assumed to be A- 
independent, as in TRG in physical phase space. Therefore all the discussion below holds 
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Figure 1: Flow equation for the self-energy in the physical space. The black circles are the full 
1-particle-irreducible vertices, while the white circles represent the TRG kernel, given by eq. (|4C|), 

only in such a case. However, the flow equation in the Bernard-Hata-Kugo approach can 
be obtained straightforwardly by including also the A-dependence of the free propagator of 
the ghost. 

Choosing as initial conditions for Z\[J] at A ^> T the full renormalized theory at 
zero temperature, the above evolution equation describes the effect of the inclusion of the 
thermal fluctuations at the momentum scale \k\ = A. At a certain scale A we are describing 
a system in which only the high-frequency modes \k\ > A are in thermal equilibrium. 
The low-frequency modes do not feel the thermal bath and behave like zero-temperature, 
quantum modes. 

We define as usual the cut-off effective action as the Legendre transform of the gener- 
ating functional of the connected Green functions Wa[J] = — Hog Z\[J]: 

, SW A [J] 
5J(-p) ' 

where we have isolated the effect of the cut-off on the free part of the effective action and 
used for the classical fields the same notation as for the quantum fields. 
The evolution equation for the cut-off effective action Ta[$] can be derived straightfor- 
wardly from eq. fl32[ ). This is done in appendix B. The flow equations for the proper 
vertices are obtained by derivation of the one for the effective action. As an illustration, 
the flow equation for the gauge boson self-energy is given in eq. (|73|) and depicted in Fig. |l[ 
We see that the recipe to obtain the TRG flow equations is very simple: one takes the 
one-loop graphs, substitutes every vertex and propagator with the full ones, apart from a 
propagator that becomes the RG kernel: 

Q 

Ka^f = «Ga iMP • At^-Goa, P t ' Ga,tv , (34) 

where the full propagator 

Ga,™ — (Goa,™ — ^A,ruj (35) 



i Tr $ ■ (A.A 1 " AT 1 ) • $ + T A [$] = W A [J] - Tr J ■ $ , $(p) = ^f^j^ , (33) 
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appears. Since the only A-dependent part in Ga, pt is proportional to the tensor defined 
in eq. fl2T|), and using the fact that is orthogonal to the remaining three tensors, the 
only component of the full propagator needed in order to compute the kernel is the one 
proportional to A^ v . It has the general form 

A w G A . pv = -A^P' A [A T (k)] , (36) 

where we defined 

Kif(k)] = Q[f(k)] + [f(k) - f*(k)] K(\k \,A) B , (37) 

and we introduced the "longitudinal" and "transverse" Feynman propagators and self- 
energies 

k z - 11l,t(k) + is 

The primes keep track of the fact that, in a generic gauge, the thermal structures appearing 
in the full propagator may correspond to an effective momentum-dependent temperature 
different from the physical one (see ref. p9j). However, multiplying the three matrices 
appearing in (|34]), it is easy to see that all the terms proportional to A^(|fc |, A) drop, and 



the only temperature appearing in the kernel is the physical one, coming from the free 
propagator. The kernel can then be written as 

K Aiflu (k) = -iA^Ao - A*)^| A8(\k\ - A) jV 6 (|fc |) B. (39) 

The above expression does not vanish in the e — > limit, and we obtain 

K A ^(k) = -A^ PT , A (k)e(k ) A5(\k\ — A) N b (\k \) B , (40) 
where e(k ) = 9(k ) — 9(—k ) and the transverse full spectral function, 

ipTA k ) = < k o) Disc A T (fc), (41) 

appears. From ([39|) we note that the product of the kernel and the inverse of the full gluon 
propagator vanishes in the e —> limit. This fact will be used in the next section to derive 
ST identities in the present formalism. 

We see that the kernel is exactly transverse and, in a covariant gauge, independent 
of the gauge-fixing parameter a. The same form holds in the TAG. It is also completely 
identical (apart from the tensor A) to the kernel of a scalar particle, which was studied in 
detail in ref. 0. 

The physical meaning of this kernel is straightforward. As the infrared cut-off A is 
lowered, the dispersion relation satisfied by the new modes coming into thermal equilibrium 
is not given by k 2 = m 2 as in the T = free theory, but is determined by the full spectral 
function. In the approximation in which one neglects the imaginary part of the full self- 
energy this simply amounts to substitute the T = mass with the thermal, A-dependent 
one. 

In section 5 we will discuss the application of the formalism outlined in this paper to 
the computation of the electric and magnetic mass in a high-temperature SU(AQ plasma. 
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4 Slavnov- Taylor identities 



In general, the implementation of the Wilson RG in the framework of Quantum Field 
Theory requires special care when gauge theories are considered. The obvious reason is 
that the introduction of a momentum cut-off on the quantum fluctuations immediately 
leads to a breaking of gauge invariance. The modified (inverse) propagator induces extra- 



terms into the ST identities, which recover their usual form only in the A — > limit JTrJ, [LI 



To see this explicitely, we follow ref. |TTJ and perform a BRS- variation (|3|) of the fields in the 
path integral of eq. (|TT| ) in order to obtain the expression generating all the ST identities. 
For definitness, we discuss the case of general covariant gauges, but the discussion goes 
closely analogously for any class of gauges. The ST identities for the proper vertices are 
obtained by derivation of the expression 



I _d*p_ I 5T A 5T A _ 5T A 5T A _ i_ < 
J (2nY \ 5Ai(v)5wU-v) Sc i (v)SvH-v) a * KP) 5c l 



A 



(27r) 4 [SAfa) Sw^-p) 6c*(p) 5v l (-p) a A 5c\-p) 

(42) 



where we have defined R A and R A (p) as the differences between the inverses of the cut-off 
and uncut-off free propagator for the gluons and the ghosts, respectively, 

K(P) = G a,V - G J > K A (p) = S A l - S" 1 , (43) 
where the color indices have been understood. 



In the first line of eq. (|42|) we recognize the usual expression for the ST identities obeyed 
by the full theory without the IR cut-off. The second line represents the breaking of BRS- 
invariance induced by the modification of the propagator. It vanishes for A = 0, as it 
should. 

Now, the main point is that while in the applications of the RG discussed in refs. [fTD|, O, 



14 ] the second line of eq. (42) gives non-zero contributions for A 0, in the present approach 
it vanishes identically, and the cut-offed Green functions obey the same ST identities as the 
ones of the full theory, the first line of eq. (|42|). This can be understood by looking at the 
difference between the R's in the two approaches. In refs. 0, [11], |14[ the RG flow describes 
the effect of integrating out quantum fluctuations (or quantum plus thermal fluctuations in 



the case of ref. [14], where the imaginary time formalism for thermal field theory has been 
considered). Therefore the cut-off is introduced by modifying the free propagator according 
to 

G, v {p) - G A ^(p) = G, u (p)e(p, A), (44) 

then 

R,Ap) = G-l(p)(Q(p,A)- 1 -l), (45) 



and analogously for the ghost propagator. As shown in ref. |TT| these expressions, inserted 



in the analogous of eq. (|42[) for the four- dimensional theory at zero temperature, give 



non- vanishing contributions such as a A-dependent gluon mass. 
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On the other hand in the present approach we are modifying only the thermal part of the 
propagators flUf) , (f25|), or (J2?D, by the substitution of eq. (|1|). Considering for definitness 
the Landshoff-Rebhan approach, we have, 

Rt = (N b (\ Po \, A) - ^(|po|))/>t,a(p) V^Cpo) Ga,^ " B • G^ 

(46) 

= (iV 6 (bo|, A) - N b (\ Po \)) Gl\ p • k A , p(7 ■ , 

where kA, po - is proportional to the kernel Ka, po -. From the discussion of the kernel in the 
previous paragraph we realize that the above expression vanishes in the e — > limit unless 
it is sandwiched between two full gluon propagators. This is not the case for eq. ( f42|) since, 
defining the amputated field-dependent functional by 

()2]U ~Wt k A (47) 



Swi(-p)5jl{p) 



we have 



5 2 Wa 

K>,ti<P h i( u j, ^ = W(bo|,A) - iV 6 (bo|)) Tr k w • G^ • WJ , (48) 

which vanishes in the e — > limit, since WJ^(p) does not contain a G l £ pfl (p) factor. The 
same is true also in the Hata-Kugo approach, where the extra contributions to the second 
line of eq. (|42D from longitudinal and gauge gluon degrees of freedom, and from the ghosts, 
all vanish by he same reason explained above. 

The Green functions defined in our approach are then gauge-invariant, in the sense that 
they satisfy the usual ST identities, given by the first line of eq. ( f42|) even in the presence 
of a cut-off A^O. 

The fact that our TRG formalism respects ST identities for every value of A is a great ad- 
vantage with respect to the other RG methods, both from the conceptual point of view, since 
a physical meaning can be given to the TRG quantities for every A (one could also think of 
possible non-equilibrium applications), and from a computational point of view, since the 
ST identities make it possible to relate propagator and vertices in a non-perturbative way. 



5 Application: the thermal masses 

The basic tool to study the behaviour of a high-temperature plasma is linear response 
theory [R0|. The perturbation of the thermal average of the vector potential induced by a 
weak source S J® is given by the well-known relation 

5Al(x) = Jdy T™ T > ab (x, y) 5J b v {y) , (49) 

where we further assume that D^dJ": = 0. 

If we consider a source that is non-vanishing in only one direction in colour space, 
then the above relation is sufficient to determine the induced perturbations on the field 
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strength operator and consequently on the gauge invariant chromo-electric and magnetic 
fields. The retarded propagator appearing in (f4"9|) is given by the difference between the 



11 and 12 components of the full propagator (|35|) , which is the object we will study. It is 
of course a gauge- dependent quantity, so the question arises of what physical information 
can be derived from it. As has been discussed in [[H], |32| the poles of the A^- and B^- 
components of the full propagator are gauge-independent and describe the exponential 
screening of chromo-electric and magnetic fields, respectively. 

In a generic gauge, the full propagator can be written as 

G A ,^ = slA^ v + hB^ + cC^ + dD^ , (50) 

where a, b, c and d are gauge- dependent 2x2 matrices, the thermal structure of which is 
also gauge-dependent [p9f] . 

Similarly, the self-energy can be generally decomposed as 

£a,^(£0 = Go^(k) - GllJk) = A^ T (k) + B^V h {k) + C^ c (k) + D^ d (k) (51) 

(the A-dependence of the S's in the right-hand side has been omitted). In general a Slavnov- 
Taylor identity will provide a relation between the four matrices in (|5"T|), so that in a generic 
gauge the self-energy will be determined by three independent functions. Nevertheless, as 
discussed in the pole structure of the temporal and longitudinal modes is identical. 



Thus there are only two independent physical modes; the spatially transverse and the 
longitudinal one. 

The longitudinal and transverse self-energies are 

E L (Ar) = B""E A ,^, 

(52) 

£ T (A;) = \A^^ V . 

The real parts of the 11 components of Sl,t coincide with the real parts of n Li T m eq. 
which will denoted as 11^ T . 

The two physical poles of the full propagator are then identified by the equations 

k 2 - nf (fc) = o , fc 2 -nf(fc) = o. (53) 

We are interested in the screening of static chromo-electric and magnetic fields. The relevant 
poles are then given by 

ml(A) = n*( ff )| , 

9=<?L 

(54) 

m|(A) = n«(g) 

<3=<?T 

where g L ,T = (?l jT = 0, |gL,T| 2 = -"Il,t)- 

The masses defined in ( |54"D are gauge-fixing-independent for any value of A, as can be 
shown by strictly following ref. [33], and considering variations in the gauge-fixing term 
which vanish for t to- 
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The evolution equations for ml and m\ can then be written down as 



A ^_ m 2 _ B ^{q)F ilv {q) 



dA i ,MW 

i+ a]?l 2 



<2=<?L 

(55) 



<?=<3T 



dA i , gn«( g ) 

where (see eqs. ( ff3|) and (|7^) ) 

F^(g) = A^-E^(g) = i y j^K A|(W (A;) [r^J^C?, -g, fc, — fc) 

-2rS pA (g, fc - g)G|' )Ar (g - k)Vf vaT {-q, k, q-k)]. (56) 

To simplify the notation, the colour structure has been understood. In appendix C we 
verify explicitly that the running does not introduce any gauge dependence. Namely, we 
show that, when contracted with B^ v (q) (A^ u (q)) and computed at q = q^ (gr), F^ v (q) 
gives a gauge-independent result. 



The system ( p5|) has to be integrated from some initial A T down to A = 0. At A we 
have to insert the zero-temperature propagators and vertices in the right-hand-sides. These 
can be reliably computed using perturbation theory, since, due to the momentum delta- 
function in the kernel, eq. ( PH) , the loop momentum is forced to be \k\ = A T ^> Aqcd 
(the external momentum is |g| = m^^^A ) <C A ). In other words, the cut-off serves as IR 
regulator for the T = theory. As we decrease A at O(gT), an IR cut-off is generated in 
the longitudinal sector (the longitudinal, or Debye, mass), whereas in the transverse sector 
a magnetic mass of 0(g 2 T) is believed to emerge non-perturbatively. 

The resummation of the thermal fluctuations from A ~ T to A ~ gT corresponds to 
the resummation of Hard Thermal Loops performed by Braaten and Pisarski 0. To go 
beyond that approximation, thermal fluctuations from gT to g 2 T have to be integrated 
out. In a perturbative framework, like HTL resummation, the different scales are not 
clearly separated, and this is reflected, for instance, in the pathologies (gauge-dependence 
and logaritmic divergence) of the next-to- leading correction to mt which require the 
ad-hoc introduction of a magnetic mass term. In the next subsections we outline how this 
non-perturbative problem can be tackled in the TRG formalism, both in covariant and 
physical gauges. 



5.1 Covariant gauges 

The class of covariant gauges is specified by taking T a = d^A a in eq. (0). The four matrices 



in (|50|) are then restricted by the ST identity: 

k"k u G A ^ = a . (57) 
In Landau gauge (a = 0) we have the stronger relation k^G^ v = 0. 
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In order to obtain an explicit form for the full propagator to be inserted in the evolution 
equations for the thermal masses, we will neglect the imaginary part of the longitudinal 
gluon self-energy, which allows us to write the full propagator as 

K,^(k) = V Pa[At] + £v Qt A L + ac 2 A 2 L A ] - aC pv Q[cA L A ] + «O p Q[A ] , (58) 

where c(k) is a gauge-dependent function. 

Due to the presence of the full three and four gluon vertices (see Fig. 1), and of the 
unknown function c, the system of coupled differential equations (|55"D is not closed. The 
knowledge of the A-dependent vertex functions would of course require the solution of their 
evolution equations which, in turn, involve higher vertex functions. 

In order to get a manageable approximation, the system has to be truncated at some 
level, in practice by using a suitable approximation of the three and four gluon vertices. A 
crucial requirement on this approximation is that it respects the Slavnov- Taylor equation 
for the full vertices to the required level of accuracy. This is essential if one wants a 
gauge-independent determination of m 2 T which are gauge-independent quantities in the 
full theory. In practice we are interested in computing the squared masses at 0{g i T 2 ). 
Looking at eq. (|55|) we see that the most dangerous a-dependent contribution comes from 
the .D^-term in the propagator of eq. (|58|). If we just use the tree- level expression for the 
trilinear vertex, then, since 

(k - p)\T t ™ x (p, -k, k-p)= g pp (k 2 - p 2 ) + p^p p - k p k p (59) 

(the group and RT indices of this vertex are trivial) , it is easy to see that this a-dependent 
term gives a contribution 0{m 2 ? ) to F /i!/ (gx), i-e of the same order as the quantity we want 
to compute. 

A way to improve the unsatisfactory approximation of tree-level vertices would be to 
enlarge the system of differential equations to include those for T^ up and T AyPvpU) setting 
to zero all the higher-order vertices appearing in the new equations. In the next subsection 
we will see that, thanks to the simplicity of the ST identities, the TAG offers a practically 
more appealing alternative. 



5.2 Temporal Axial Gauge 



The ST identities for the three and four gluon vertices in the TAG are 



^a6c 



(p, q, r) 



f abC { Gl](q) S 13 - 



lt£«r(P,ff,r,a) 



ljki 



8 11 



(60) 



r i HI r i lkl 

+*f ace ri 6 i> + r, q, sj\ 5 lk + if ahe Tl cd (p + q, r, S ) . (61) 



In (p0| ) we have defined G A ^ u (q) = G A ^ pu (q) — u p u u /a Q[l], which is well-behaved in the 
TAG limit a 0. 
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Neglecting again the imaginary part of the longitudinal self-energy, the full propagator 
can be written as 

G A ,^ = -V Pa[At(*0] " Q' [H^A h (k)} , (62) 
where H^ v has been defined in (PU|). 

The self-energy is given by 

^u{k) = V P A 1 [A T (A;)- 1 - k 2 ] + B, v Q' 1 fA L (A;)- 1 - k 2 } . (63) 



Now, as promised in subsect. 3.2, we show that, as a consequence of the ST identity 
the spurious poles in l/k' and l/(k' ) 2 appearing in H^ v do not give any contribution to the 
running of the two physical masses m\ and m\. Indeed, from ([$(]) we see that these poles 
appear as coefficients of the tensors C^k) or D IMl ,(k), and then they are always multiplied 
by k^ or k v . So the analogous singularities appearing in the propagator G^\ T (q — k) in 
eq. (|56|) are proportional to (q — k)\ or (q — k) T . Applying the ST identity (|60D, at least one 
of the two full three-gluon vertices in Q55D can be written as the difference of the inverse 
propagator G _1 at the loop momentum k and the one at the external momentum q. The 
first one vanishes when multiplied by the kernel K.\ piJ (k). The second one vanishes at 
q = (ft, if F^ u (q) is contracted with B^ u (q) and at q = q T if it is contracted with \A^ V (q), 
so it does not contribute to the evolution equations for m£ and m\. The cancellation of all 
spurious singularities can also be proved for the denominators in (B^) in a closely analogous 
way. 



In the past, the simple form of the ST identities (p0|) and (|61|) has motivated the choice 
of the TAG in non-perturbative studies of the infrared limit of the gluon propagator, both 
at T = |34| and at T ^ |35], [36|. These studies were based on truncations of the exact 



Schwinger-Dyson equations for the self-energy. In |34] and [Bg] the observation was made 



that the transverse part of the full three-gluon vertex, defined by 

-i(T)o6c 



PAri T il c (p,g,r) = o, (64) 



and the requirement that be free of kinematical singularities, vanishes when any 

of the external momenta goes to zero. Then, in the infrared regime of interest here, the 
relevant part of the three-gluon vertex is given by the longitudinal one, which is completely 
determined in terms of the inverse full gluon propagator by means of eq. (|60D \ $% . By 
expressing the four gluon vertex in terms of the longitudinal three-gluon one, the set of 
infinitely many Schwinger-Dyson equations is turned into two coupled equations for 
and m\. 

In the present framework the same approximation of the A-dependent three and four 
gluon vertex can be implemented into the exact (non-truncated) equations (|5|). Since the 
vertices are fully determined by Hh{k) and Hr(fc), the possibility arises of closing the system 
of differential equations for the thermal masses. The simplest approximation would be to 
take n L)T (/c) ~ mL T , in which case one is left with only two coupled ordinary differential 
equations for m\ and m\. 

This is a good starting point for a non-perturbative computation of the screening masses. 
In a forthcoming paper we will discuss this and more refined approximations and present 
the results of a numerical computation. 
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6 Conclusions 



In this paper we have discussed the extension to gauge theories of the Wilson RG-inspired 
formalism introduced in ref. ||. The main attractive feature of this formalism is given 
by the fact that introducing the infra-red cut-off on the thermal sector of the theory does 
not break gauge invar iance. Namely, even in the presence of a non-zero thermal IR cut-off 
A, the ST identities are formally the same as at T = 0, provided the Green functions are 
interpreted as thermal averages with respect to a density matrix p(A). 

This provides the basis for a gauge-invariant renormalization flow, in which the running 
does not introduce extra gauge dependence. In computing a gauge dependent quantity, the 
result at the end of the running (A = 0) will be the finite temperature counterpart of the 
same quantity computed at T = in the same gauge. On the other hand, the TRG running 
for a properly defined gauge independent quantity will be completely gauge independent. 

From a practical point of view, a gauge invariant RG flow presents the great advantage 
that ST identities can be used as a guide to reduce the number of independent invari- 
ants to be included in approximation schemes based on truncations of the exact evolution 
equations. 

The need for a proper control of gauge invariance in computation methods for the 
thermal masses in hot QCD has recently been emphasized by many authors (see for instance 
f38fl). We have discussed this problem in our formalism, pointing out that an approximation 
of the full set of TRG equations based on lowest-order truncations, in which only the 
running of the two-point functions is included, leads to a gauge-dependent effect of the 
same order as the subleading corrections we are interested in. Due to the simplicity of the 
ST identities, the TAG turns out to be a more appealing gauge, in which improved vertices 
can be defined in terms of inverse propagators and a closed set of approximate equations 
can be obtained. Remarkably, the HTL propagators and vertices are gauge independent 
and satisfy the same ST identities as those of the TAG |y. Then, in stepping beyond 
HTL rsummation, one could impose that the same identities are satisfied also by the new 
propagator and vertices, even in covariant gauges. 

Work along these lines is now in progress and will be presented in a forthcoming paper. 
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A TRG generalization of Hata-Kugo proof 

The partition function flTBD can be rewritten as 

Z{0) = TrfPe-^p] , (65) 
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since P projects onto the subspace with FP number 0. Then we notice [26 that P can be 
written as 

P = 1 - {Qbrs, R} , 

where R is a suitable operator of FP number —1. By inserting this in fl65|) we find that the 
contribution containing R vanishes 

Tr[{Q B R S , R}e~^p] 

= Tr[R e~ wQc p Qbrs) + Tr [R Qbrs e~~ nQc p] (cychcity) 

= Tr[Re-^ c Q B Rsp] +Tr[RQ B R S e-" Qc p] (since [Q B rs, p] = 0) 

= Tr[R {e~* Q % Q BRS } p) = (since [Q BRS , Q c ) = iQ BRS ) 

The crucial point is that the non-equilibrium density matrix p commutes with the BRS 
charge, due to the adiabatic switching-off of the interaction at t — to. 



Therefore one finds (14). 



B TRG equation for the cut-off effective action 



In this appendix we derive the TRG flow equations for a general gauge. 



By deriving eq. (|33|) with respect to A and using (|32|) we get 



(66) 



In the RHS of this equation there is a field-independent piece coming from the full propa- 
gator. We isolate it and define a new functional f : 



(27T) 



J 2 W A [J] 

SjJj v 



(27r) 4 Ga,^ - G 



A,fMp ' I^A, per ' G A ,<tu 



(67) 



where A is the full propagator (|35f) . In this way the evolution equation ( pq) becomes 
(dropping the field-independent term) 







A ^X r A[ $ ] = ^Tr{K A ^-f A ^} 
where the kernel K is given in (34). 



(6? 



We now have to express the functional in the RHS of flog) in terms of From the 
definition (03) one has 



- J M = Gqa" ■ A v + (2vr 



6A V 



Deriving this relation with respect to j„, it follows 



(27T)* 



5 2 W A 



G OA,lp + ( 27r ) 



t ^ 2 r A 

8A v 8Ar 



9w 



(69) 



(70) 
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We now isolate the self-energy 



(2 " )8 ^^ = (27r)4SA - + r "- [$] - 



Substituting eqs. (|6"?D and ( |7T|) into (0), one finds 



rint 
/ 



(71) 



(72) 



The evolution equations for the various vertices can be found by expanding in powers of 
$. As an example we give the flow equation for the gauge boson self-energy (see Fig. [I]). 
By deriving QBBD with respect to A*(p), A^(— p), and setting the fields to zero, we get (in 
the following the gauge group indices a, b, c, d are explicit, and no summation is made on 
the index a) : 



A 



dA 



(p) 



ii 



(2tt; 



4 2 A,per 



(*) 



i A,p CT 



^- i (p)M^ 1 (-p) 



where the vertex of T is found from 



A,p<r 



<5^' 1 (p)5X' i (-p) 



a,l, 



paa66,llij'/ j j \ 

1 A,^pcr IP? — Pi ~ K i K ) 



$=0 



-2 rJ2j?(p, *: - p) [GA,Ar] fc/ (p - A;)r A b ^(-p, k,p — k) 
where we have the three and four gluon vertices 

5 3 r A [$] 



ifcZ 



^abc,ijk 
A, pup 



(p, q, r) 



-r^abcdAjkl / \ 

T A,aupi (P,q,r,s) 



8A*{r)8Al{q)8Afa) 

5 4 r A [$] 



5A l a (s)5A k (r)5Ai(q)8Ai(p) 



(73) 



(74) 



C Gauge-fixing independence of the TRG equations 

for ml and m\ 



We are interested in the effect of a change of gauge-fixing function: 

on the evolution equation for the self-energy. 

Following ref. [33 we will start from the basic relation 

AW JT [J]=-J i A Xi , 



(75) 



(76) 



where AW^J] is the change in the generating functional W A [J] under the transformation 
Q7§), and 

A X i = (D l -M- 1 -AF). (77) 
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Under the same change (|75|) the evolution equation for the effective action T A [$\ is 
modified by the extra terms 

A^Ar A [$]=#i[$] (78) 



(79) 



where 

Him = -WlvjjA^AG^.^}. (80) 

Deriving the above expression twice with respect to and A u , we get the variation in 
the TRG running of the self-energy £a,^(<?) induced by ([751). Using (|76|) , and after some 
algebra, we find that the first contribution ([T^]) can be rewritten in the form 

^(27r) 8 Tr {G 0A \ ia (q) ■ K A>ij (k) ■ X ijav (k, -k, q, -q)} + (// <- u) . (81) 

To obtain the evolution equations for Sl(q') and Sx(g), we have to contract with B^ v (q) 
and ^A^q) respectively (see eq. (|52"|)). The inverse full propagator at go = in a generic 
gauge can be written as 

G-l] 11 (g = 0, q) = ^At 1 + B m A L 1 + Dfjtafd ■ (82) 

The structure is absent at go = due to T-invariance, and moreover 

ImS(g = 0,g-) = 0. 

The contribution of ([FjJ) to the variations in the running of S L X and S" * s then given 

by 

A^-AE^go = 0,q) = z(2tt) 8 A L \q)B au {q) Tr{K Aiij ■ X ijau } , 

d (2ir) 8 
A— AS" (go = 0, q) = i^LA T \q)A au {q) Tr {K Aij ■ X ljau } . (83) 

From the above expressions we see that if we compute El and St on the longitudinal and 
transverse mass shells, respectively, the gauge dependence of (^) vanishes. This conclusion 
could be spoiled by the presence of a pole of the function Xij ai/ for q on the same mass shell. 
As discussed in |33j the pole structure of this function is determined by the poles of the full 
ghost propagator, which are gauge-dependent. As a consequence this coincidence between 
the gluon and the ghost poles could take place only accidentally, in some particular gauges. 
If only one gauge in which this is not the case exists, then the gauge independence of the 
gluon dispersion relation can be ascertained. 

The second contribution to ([T8|), eq. (|8"0|), is absent in the Landshoff-Rebhan formulation 
of covariant gauges as well as in the TAG, so the proof is complete in these two cases. We 
will then consider this contribution in the only case of practical interest in which this is 
different from zero, that is in covariant gauges a la Hata-Kugo. The only gauge- dependent 
part in the tree- level propagator is proportional to the tensor D p<T {k): 

A |v AG ^ = ~ 2l£ ^t^ D ^ k)B ■ (84) 
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The contribution (]80f) then gives 

(2n) 8 e^^Tr {D pa (k)G A>pp ,(k) • B • G w (fc) • fS^fe, "9, *, "*)} , (85) 

where the vertex T°£ b pupl(T , is the one defined in (j74|). 
Now we use the ST identity to write 

- ik p G A , p(T (k) = aT { ${k) ■ B' ■ D A (fc) , (86) 

where T) A (k) is the full ghost propagator, 1^ ^ the w — c vertex, w a being the BRS-source 
defined in eq. ([?]), and B' = diag(l, — 1). Contracting this vertex with the four-point 
function in (|85|) and using the ST identities for the three and four gluon vertices, the 
parenthesis in (|85| ) turns out to have the structure 

G^ a (q)Z au (q) + Gl] va {q)ZaM) + G^ a (q) G^ x (q)Y aX (q) , (87) 
where G^^q) = G A ^ Q (g) - -A; M A; Q Q[1]. 

Contracting again with B pv (q) and with iy4 M ^(g), we conclude that also the contribution 
of ( |30D vanishes when q is on the (A-dependent) mass shell. 
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